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ABSTRACT 

We develop a novel method to separate the components of a diffuse emission process 
based on an association with the energy spectra. Most of the existing methods use some 
information about the spatial distribution of components, e.g., closeness to an external 
template, independence of components etc., in order to separate them. In this paper 
we propose a method where one puts conditions on the spectra only. The advantages 
of our method are: 1) it is internal: the maps of the components are constructed as 
combinations of data in different energy bins, 2) the components may be correlated 
among each other, 3) the method is semi-blind: in many cases, it is sufficient to assume 
a functional form of the spectra and determine the parameters from a maximization 
of a likelihood function. As an example, we derive the CMB map and the foreground 
maps for seven yeas of WMAP data. In an Appendix, we present a generalization of 
the method, where one can also add a number of external templates. 
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1. Introduction 

In many observations, understanding the diffuse emission components plays a significant role in 
the interpretation of results. For instance, diffuse background modeling is crucial for the observation 
of the CMB fiuctuations in radio frequencies and for constraining DM annihilation in gamma-rays. 

Diffuse emission in a pixel i at an energy can be described as a sum over emission components 

m 

where 5^ is the line of sight density of component /x at pixel u]^ is the corresponding energy 
spectrum at energy (the matrices u are also called mixing matrices), denotes random noise 
(instrumental or physical). 

In various applications it is important to find the energy spectrum and/or the spatial distribu- 
tion of the components. In many cases, one can neglect the dependence of the spectrum on pixels 
within the region of interest, i.e. u]^ = u^. Let us also for the moment neglect the noise term r^. 
Then the problem is to find and 5^ from the system of equations 

m 

11=1 

Suppose that the total number of energy bins (frequency bands) k is larger than the total number 
of significant components of the emission m. We also assume that the number of pixels is larger 
than both k and m. Then, in general, this system is both overdetermined and redundant. It is 
overdetermined because the rank of the matrix rf^ is k while the rank of the matrices and 5^ is 
m. It is redundant because there is a symmetry of multiplying the matrix u from the right with a 
matrix A and multiplying the matrix s from the left with the inverse of A 

u — > u ' A^ s — > A~^ • s. (3) 

Different methods can be classified according to the usage of the overdetermination of the system 
to most efficiently extract the components of the emission and how they deal with the redundancy 
of the system. 

One of the most simple cases is when either the spectra u or the spatial distributions s of the 
components are known. If u is known then s can be found from Equation ([2]), and vice versa, 

s = {u^u)~^u^ 
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Denote / = (u^u) ^u^ . The matrix / satisfies / • 2i = /, it is the Moore-Penrose pseudo-inverse 
matrix for the rectangular matrix u. Analogously, h = s^{ss^)~^ is the pseudo-inverse for 5, s-h = I. 

If the uncertainty does not depend on both the energy and the position, = a, then Equations 
([4]) can be also derived by minimizing 



X 



(5) 



There is a summation over the repeated index ji. In the following we will assume a summation over 
the repeated indices unless otherwise stated. 

In general, the methods to separate the diflFuse emission components have the following steps: 

1. Choose a data analysis principle. 

2. Find a functional such that the principle from p.l is realized at the minimum (maximum). 

3. Break the redundancy which is not broken by the functional. 

In many methods the separating principle is the maximum likelihood. In the case of Gaussian 
statistics, the minimum of the is the maximum likelihood point. The has the change of the 
basis degeneracy (Equation ([9])). Some examples of the ways to break the degeneracy are as follows 



(a review of the methods can be found in Leach et al. (2008)) 



1. Independent components analysis (ICA), Hyvarinen (1999); Maino et al.| ( |2002| ). There 
is an additional term in the functional that depends on the mutual information (the minimum 
corresponds to maximal independence). 



2. Maximal entropy method (MEM), [Hmshaw et al.| ( |2007D ; [Hobson et al.| ( [l998[ ); [Stolyarov 



et al. (2002). There is an additional term that quantifies a closeness to some prior (usually 



an external template). 



3. Generalized morphological component analysis (GMCA), |Bobin et al. (2007). There is 
an additional term that quantifies the separation of components in a certain dictionary. 



Other methods have a diflFerent physical principle to start with. For example. 



1. Principal components analysis (PGA), e.g., Francis & Wills (1999); Steiner et al. (2009). 



The functional is the variance in the space of data vectors. The first PC A component corre- 
sponds to the direction of the largest variance, the second component corresponds to the next 
largest variance direction etc. The degeneracy is broken by choosing the components in the 
order of decreasing variance. 
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2. Internal linear combination (ILC), [Eriksen et"aL (2004); Hinshaw et al. (2007). One uses 
the fact that the expected variance of a sum of two components is larger than the variance of 
either of the components. The CMB map is modeled as a linear combination of data vectors 
at different WMAP frequencies with the constraint on the coefficients ensuring that the CMB 
signal is preserved in the linear combination. The CMB map is determined by minimizing the 
variance across the sky of the linear combination map. 

3. Correlated component analysis (CCA), [Bedini et aL (2005). The basic observables are 
not the fluxes but the two-point correlators. One uses some parameterization of the source 
correlator matrices and the energy spectra of the sources. The parameters are estimated from 
the least squared difference between the observed correlators and the model correlators. The 
redundancy is broken by choosing a particular functional form of the spectra and/or assuming 
that certain components are uncorrelated. 



On general grounds, the degeneracy can be broken by either putting constraints on the energy 
spectra or on the spatial distribution of the components. Most of the models (apart from the CCA) 
assume some properties of the component distributions. Usually this means either using external 
templates or assuming some independence of the components. In CCA one uses the functional 
form of the spectra to break the degeneracy. This allows one, in principle, to separate correlated 
components internally. 

In this paper, we describe a novel method, which we call the spectral components analysis 
(SCA). In this method, we use the maximum likelihood for the fluxes and assume a functional 
form of the energy spectra to break the redundancy. In comparison with the methods that put 
constraints on the spatial distribution of the components, the new method allows one to study 
correlated components using only internal data. In contrast with the CCA method, one formulates 
a model directly for the fluxes rather than for the correlators: the decomposition is simpler and the 
calculations are easier in this case. 



2. Spectral components analysis 

In this section we formulate the SCA method to separate the components of a diffuse emission 
process. The components of the emission are modeled by linear combinations of the data in different 
energy bins. The degeneracy between the components is broken by assuming a functional form of 
the energy spectra. The linear combination parameters and the parameters of the spectra are found 
from a minimization. 
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To start, we introduce vectors with components labeled by pixels, e.g., the data in a bin with 
the central energy is written as rfc^ = {rf^} where index i labels the pixels. In vector notations. 
Equation ([T]) takes the form 

da = ^K^fi + ^a, (6) 

where we assume that the energy spectra do not depend on pixels. We assume that the num- 
ber of the emission components is smaller than the number of energy bins and model the spatial 
distribution of the emission components as linear combinations of data vectors 



s 



Following the WMAP terminology (Hinshaw et al. 2007), we will sometimes call the emission 
components determined as linear combinations of the data vectors by internal linear combination 
(ILC) vectors. The spectra u and the linear decomposition coefficients / are found by minimizing 
the 

^^yiAz<mi. (8) 

— ^ 

The matrices u and / can be determined only up to the change of basis degeneracy, which corre- 
sponds to multiplying u from the right by a matrix A and multiplying / from the left by the inverse 
matrix 

u-f = uA-A-'f. (9) 

We break the degeneracy by assuming a functional form of the spectra parameterized by a set of 
parameters q 

u''^{q) = u^{E^;q). (10) 

Before we formulate a general algorithm, let us consider the simplest case. We fix the parameters 
of the spectra q and assume homogeneous independent of energy uncertainty, a. Then according to 
the first equation in Q 

/ = {u^u)-^u^. (11) 

The product of the linear combinations matrix defined by this equation and the matrix of spectra 
equals the unit matrix 

E/x = '^M' (12) 

a 

which means that a component with spectrum contributes only to the linear combination and 
no other components. In general, due to different properties of noise in different energy bins, the 
conditions in Equation (12) will be slightly violated. 
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Notice, that depends quadratically on / and non-linearly on q 



x\fA) = Y.— 2 ' ^ ■ (13) 



a' 



The algorithm has two steps: 

• Find / from quadratic fitting. 

• Minimize with respect to g by a non-linear fitting. 

At each step of the non-linear fitting procedure we choose a set of parameters q that describe 
the spectra and find the best fit / by differentiating with respect to / and by solving the 
corresponding linear in / equations 

• d^)9J^ = ^<{da • rfik, (14) 

a a 

where by (rf/^ • dj)g^ = gaijd'^pdi^ we denote the scalar product in the space with metric g^. In our 
case, it is the inverse of the standard deviation squared 

gaij = ~~2^ij- (15) 
^ a 

In general, gaij is the inverse of the covariance matrix. 



Equation (14) is a usual linear equation, the non-trivial part is that the unknowns /J have two 
indices instead of one. We can solve this equation, if we choose a single index corresponding to all 
possible pairs of indices. The mathematical structure behind this operation is the tensor product, 
i.e., we represent the coefficients as vectors in a linear space F which is a tensor product of two 
linear spaces: the space D spanned by vectors df^ and the space S spanned by vectors Sj^^ F = D<S)S. 
Suppose that the dimensions of spaces D and S are k and m respectively, then the dimension of F 
is k • m. The basis vectors in F are denoted as 

U = d^® s^. (16) 

Loosely speaking the operation of the tensor product is a substitution of a pair of indices /x} by 
a single index k. Denote k = and A = {7, z/}, then the coefficients /J can be represented as 



a vector with a single index /a = /J- The right hand side of Equation (14) is a vector 

v^ = J2<A-d^)9.- (17) 

a 

The sum on the left hand side is a matrix 

M'^' = J2«idrd,)g^- (18) 
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In the tensor product notations, Equation ([14|) takes the form 

^M^^fx = V\ (19) 

A 

The solution is 

fx = Y,M-\,V\ (20) 

hi 

This equation gives the minimal value of foi" ^ given set of parameters q. If we substitute 
fj{q) = fx{q) in Equation (13), we get x^{q) a function of q only. Now we can use a non-linear 



minimization procedure to find the best fit parameters q^. 

We note that the spectral components analysis allows one to find both the spectra and the 
spatial distributions of components with only two rather mild assumptions: 

• The spatial distribution of components does not depend on the energy. The only dependence 
comes in the form of overall normalizations, i.e., energy spectra. 

• These energy spectra have some functional form. 



3. Example: WMAP data 



In this section we use the SCA method described in the previous section to study the WMAP 
data ( Jarosik et al.|[20TT ). The WMAP experiment measures the fiuctuations of temperature in five 
frequency bands with the center frequencies (in GHz) 



K Ka Q V W 
22.8 33.0 40.7 60.8 93.5 



(21) 



The data are naturally split into five vectors corresponding to the five frequency bands. In the 
analysis we use seven years of data (Jarosik et al. 2011) with pixel size corresponding to the 
HEALPix parameter nside = 512 (Gorski et al. 2005). We average the data to get ^ 1° pixel size 
resolution (nside = 64). We mask the Galactic plane within ±30°, detected point sources (Gold 



et al. 2011J), and the one hundred brightest pixels in the K band which we consider to be outliers. 
The remaining number of pixels is 24,228. 

In the analysis, instead of using the data itself in the formulation of the models in Equation 
([t]) we use the data vectors randomized with the detector noise. In this case, the expected does 
not depend on the number of components. Without randomization, the expected decreases as 
the number of components increases. For example, in the extreme case of a five-component model 
for the five frequency bands the x^ is zero. 
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We start the analysis by assuming only one component of the emission. The index of the energy 
spectrum in this case is n = —0.13. The reduced is /dof ^ 14, which is relatively large. We 
conclude that one component is not sufficient to describe the WMAP data. 

For two components, we get the following indices of the energy spectra rii = 0.04 and 712 = —3.1. 
The fit in this case is sufficiently good, /dof ^ 1.16. Based on the spectra, we conclude that 
most of the contribution to these components comes from the CMB and the synchrotron emission 
respectively. These are the two strongest components at high latitudes. 

Note, that in this approach we don't need to assume that the CMB signal exists. Instead, we 
derive from the data that there is a component with an index n ^ which can be interpreted as 
the CMB ffuctuations. 

The two-component model is already very good. If we add a third component and allow the 
indices to be unconstrained, then the fit converges to some collinear vectors. In this case, one may 
choose to use some of the external knowledge. In this analysis, we assume that the third significant 
component is the dust emission and fix the index of the third component n^ust = 2, which is a usual 



assumption for the spectrum of the dust emission (e.g., |Hinshaw et al. 2007). Although it is true 



only approximately ( Finkbeiner et al.||1999 ), the exact value is not important since this component 



is subdominant to the CMB and the synchrotron emissions at the WMAP frequencies. We also fix 
the index of the CMB emission ncMB = 0. The reduced is x^/dof ^ 0.97. The best-fit index of 
the second component is 712 = —2.6. We interpret this component as a combination of synchrotron 
and free-free emission. 

The components are modeled by taking linear combinations of the WMAP temperature maps. 
The best fit linear combinations are: 

/cMB = (-0.37 0.46 0.67 0.34 -0.11 ) 

/^yn+ff = ( 0.84 -0.18 -0.42 -0.28 0.04 ) (22) 
/dust = ( 0.31 -0.44 -0.52 -0.01 0.66 ) 

Note, that the sum of the components for the CMB vector is approximately equal to 1 

E/cMB = 0.99. (23) 

(3 



In the ILC method (Hinshaw et al. 2007) this is a constraint on the linear combinations describing the 
CMB map. In our approach this is a property of the linear combination of data vectors corresponding 
to a component with n = spectrum. The sum of components for the foreground models is less 
than 0.01, i.e., the CMB signal is canceled in these models. 
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The WMAP temperature maps and the residuals after the three-component model subtraction 
are shown in Figures [T] and [2} The residuals are sufficiently small with little large scale structures. 
The SCA model of the CMB map is presented in Figure [3| It is also compared to the WMAP 



internal linear combination (ILC) model of the CMB map ( [Gold et al. 2011). In Figure [4] we 
present the map of the component corresponding to the spectral index n = —2.6. We compare it to 
the sum of the synchrotron and the free-free emission WMAP Markov chain Monte Carlo (MCMC) 



models ( [Gold et al.1|2011p . In Figure [5] we compare the SCA dust model to the WMAP MCMC 



dust model (Gold et al. 2011). Although the SCA models and the WMAP models of the emission 



components are derived using completely different methods, there is a fairly good agreement among 
the models. 

In Figure [6] we compare the angular power spectra of the SCA models to the corresponding 
angular power spectra of the WMAP MCMC models. The CMB angular power spectra are very 
similar to each other, while the spectra for the SCA models of the Galactic emission have more 
power at small £ relative to the WMAP MCMC models. In Figure [7| we present the root mean 
squared (RMS) of the temperature maps for the three SCA components at the WMAP frequencies. 

The SCA models generally reproduce the WMAP emission models. The SCA models have 
a higher level of random noise, which is expected for the models derived from internal data in 
comparison with the models that also use external data. The only assumptions that were used to 
derive the SCA maps of the components are the power-law energy spectra with fixed indices for the 
CMB and the dust emission, while the index for the third component is derived by minimizing the 



4. Conclusions 



In the paper we propose a novel method to separate the components of a diffuse emission 
process based on an association with energy spectra. The new method is 

• internal (the components are modeled as combinations of data at different energies); 

• semi-blind (we assume the functional form of the spectra and fit for the parameters); 

• can separate correlated components. 

Using only internal data we avoid systematic uncertainties related to the usage of external data. For 
instance, in various diffuse emission models one needs to know the distributions of atomic, molecular, 
and ionized gas. The distribution of the atomic hydrogen HI (Hi) is traced by the hyperfine 
splitting emission at 21 cm (Kalberla et al. 2005; Kalberla fc Kerp] 2009). The corresponding 
observations have uncertainties related to the determination of the radial velocity of the gas and the 



spin temperature Ts ( [Dickey et al.||2009 ; Johannesson et aL]|2010 ). The distribution of the molecular 
H2 gas is traced by the emission from the CO molecules ( Dame et al.||2001 ), some uncertainties may 
arise from a space dependent conversion coefficient Xco (e.g., Strong et al.||2004 ). The distribution 



of ionized hydrogen HII can be modeled by studying the H- alpha emission (e.g., Finkbeiner 2003 



Gaensler et al. 2008). In addition to uncertainties in the distribution of the interstellar gas, the 



models of cosmic ray (CR) propagation in the Galaxy (Ginzburg & Syrovatskii 1964; Strong et al. 



2007; Vladimirov et al. 112011; Cholis et al. 2011) have uncertainties in the distribution of the CR 



sources, in the propagation parameters, and in the distribution of the interstellar radiation (e.g.. 



Moskalenko et al. 2007| ). In general, the uncertainties arise from the fact that one either needs to 
make some assumptions when using a certain map as a template for an emission process or there is 
some processing involved in constructing an emission component from external data. 

The advantage of internal methods is based on the following simple observation: if a physical 
process contributes significantly to the data, then one should be able to use the data to describe the 
spatial distribution and the energy spectrum of this emission process. Otherwise, if the process is 
insignificant for a particular observation, then one does not need to include it in the consideration. 

In the SCA method, one does not need to assume a priori either the number of the emission 
components, their spectra, or the spatial distributions. All this information is obtained by fitting 
the model to the data. As a result, the SCA method is a unique tool to search for new components of 
emission with spectra sufficiently different from the astrophysical foregrounds. In template fitting 
and in galactic propagation frameworks, a new component can be observed as a residual. The 
spectrum and the significance of a "residual" component are usually biased. An advantage of 
the SCA method is that the new component and the galactic components are fitted to the data 
simultaneously, i.e. all components in this approach are on equal footing and the there is no reason 
to expect a bias for any of the components. 

Another advantage of the SCA method is the possibility to separate correlated components. 
This property distinguishes the SCA method from other internal data analysis methods such as 
internal linear combinations, principal components analysis, or independent components analysis. 
One of the disadvantages of the SCA method is that we need to assume a homogeneous scaling of 
the components, i.e., that the spectra do not depend on the position in the sky. Although it is not 
true in general, it is usually a good first order approximation, especially if the region of interest is 
small. 

In the paper, we use the SCA method to study seven years of the WMAP data. We find that 
at high latitudes a three component model describes the data with a sufficiently good accuracy. We 
interpret the three components as the CMB, the thermal dust, and a combination of the synchrotron 
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and the free-free emission. The derived spatial distributions of the components are generaUy con- 
sistent with the corresponding spatial distributions provided by the WMAP collaboration. 

The SCA method is sufficiently universal. It can be applied to other types of diffuse emission 
data, such as infrared, optical, x-ray, and gamma-ray. One can also use any representation allowing 
linear space interpretation, e.g., coordinates space, spherical harmonics, wavelets. The method has 
a straightforward generalization to the case where some of the components are modeled by external 
templates (Appendix [A]) . 
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A. General case 

In this appendix we present an algorithm that includes some number of external templates 
together with components modeled as internal linear combinations of data vectors (ILG vectors). 
In the most general form, we will not put any constraints on the spectra for some of the ILG vectors. 
As we have discussed in the introduction, if there is more than one ILG vector for which we do not 
put constraints on the spectra, then one needs to break the change of the basis degeneracy (Equation 
([9])). In the algorithm described below we break this degeneracy by choosing some sub-matrix of 
the linear combinations matrix / (Equation ([7|) to be the unit matrix. Gomponents for which we 
do not put constraints on the spectra have an independent scale factor in every energy bin. The 
process of subtracting these components is called marginalization. 

In the general SGA algorithm we can consider the following models for the components together 
with the notations (bold font letters denote vectors, while usual font letters denote numbers): 

1. Gomponents with known spatial distributions (templates) that we marginalize over, t5^e^ (no 
summation) . 

2. Gomponents with known spatial distribution for which we assume a parametric form of the 
spectra, w{p)^ejy (no summation). 

3. Gomponents with unknown spatial distributions (modeled as ILG vectors) that we marginalize 
over, u^f^V/s (no summation over /i). 
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4. Components with unknown spatial distributions (modeled as ILC vectors) for which we assume 
a parametric form of the spectra, u^{q)f^^[3 (no summation over /x). 

The general has the form 



W - <(g)/^/3 - </^V/? - w{p)le, - <e. , (Al) 
— I gcx 



where |v|^^ = (v,v)^^ is the norm in space with metric Qa- The metric is given by the inverse of 
the standard deviation squared, or, in general, by the inverse of the covariance matrix. 

The algorithm has the following steps: 

1. Marginalize over This is equivalent to projecting the vectors dc,, and onto the space 
perpendicular to the space spanned by e^. The residual has the form 



' 9a. 



2. Choose parameters p. Subtract w{py^e^ from and v^. 

3. Choose and marginalize over the parameters by projecting d^ on /^v^. 



4. Choose parameters q and marginalize over (Equations (13) to (20)). 



5. Repeat steps 2, 3, 4, (and 5) to find the best fit non-linear parameters p*, g*, and 

Given the best fit parameters p*, g*, and the maps of the components are reconstructed as 
follows: 

1. Components with templates and constrained energy spectra are given by t(;(]9*)^e^ (no sum- 
mation). In order to get the other maps, we subtract these components from the data and 
define d'^ = d^, - w{p^y^e^. 



2. The maps for the templates with marginal spectra are obtained by marginalization of d^ with 
respect to the templates e^. The scaling coefficients are obtained from minimizing the 



< = EK'e,k^-^"^ (A3) 

where E~^^^ is the inverse of the matrix of scalar products E^y = {e^^ey)g^. The maps are 
given by w^^^y (no summation). 



3. The ILC components are obtained by linear combinations of model vectors derived from 
the data vectors d'^. If Vc, = dc,, then = v^, — w{p^)^ejy. We also marginalize over the 
templates with unconstrained spectra by projecting d'^ and on e^: d'^ d'^^ and 
v^^. The maps for the ILC components with marginal spectra are given by u^J^fj^^'^ 
(no summation over /x). The overall coefficients are found similarly to w^^ by minimizing 
the residual that depends on d'^^ 



(A4) 



where = /f^v'^ are the ILC combinations and F ^'^^ is the inverse of the matrix of scalar 
products F^^ = {if,X)gc.- 

4. The maps for the ILC components with functional forms of the spectra are equal to u^{q^) f^^^ 



(no summation over /x). The coefficients are found from Equations (14) to (20), where we 
substitute ^^{q^) dc, d'^. 



The general SCA method gives simultaneously the spectra and the maps of all four types of models 
for the emission components: templates with and without assuming a functional form of the spectra 
and ILC combinations with and without assuming a functional form of the spectra. 

We note that in some cases instead of a non-linear fitting procedure one can use a convergent 



iterative process (e.g. Bobin et al. 2007; Tsalmantza & Hogg 2012) to define the spectra and the 



spatial distributions of the emission components. 




Fig. 1. — Left: maps of seven years of WMAP data in K, Ka, and Q frequency bands ( Jarosik et al.||2011 ). 
Right: residuals after subtracting the SCA models (Figures |3| [ij and [5]) divided by the instrumental noise. 




Fig. 2. — Same as in Figure [T] for V and W frequency bands. 
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Fig. 3.^ Top: the SCA model of the CMB fluctuations. Middle: the WMAP ILC model of the CMB 
( [Gold et aT1|201lD . Bottom: the difference between the SCA model and the WMAP ILC model of the 
CMB divided by the instrumental noise. The SCA CMB model is very similar to the WMAP CMB model, 
but has a little higher random noise level. 




Fig. 4. — Top: the SCA model of the synchrotron+free-free emission. Middle: the WMAP MCMC model 
of the synchrotron+free-free emission (Gold et al. 2011). Bottom: the difference between the SCA model 
and the WMAP MCMC model divided by the instrumental noise. 




Fig. 5.— Top: the SCA model of the thermal dust emission. Middle: the WMAP MCMC model of the 
thermal dust emission ( Gold et "aL||2011 ). Bottom: the difference between the SCA model and the WMAP 
MCMC model divided by the instrumental noise. 
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Fig. 6. — Angular power spectra for the SCA models and the WMAP models presented in Figures |3| |4j 
and [5j The angular power spectra of the CMB are very similar to each other. The SCA models of the 
synchrotron+free-free and the dust emission have a little more power at small I than the WMAP MCMC 
models. Large up and down fluctuations of the foreground emission Q's is due to the symmetry properties 
of the mask. In calculating the Q's, we divide by the fraction of the sky that is unmasked. 
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Fig. 7. — Root mean squared of the temperature maps for the SCA models (Figures |3l [4| and [5]) at 
different frequencies. We assume power-law energy spectra. The indices for the CMB and the thermal dust 
emission are fixed, ncMB — and n^ust — 2. The index of the third component and the normalizations 
are found from fitting the models to the data in Section [3j The index for the third components is found 
to be ris+ff = —2.6. This component is interpreted as a combination of the synchrotron and the free- free 
emissions. 
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